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ABSTRACT 

We generalize a formula for the Goldreich- Julian charge density (p G j), originally de- 
rived for rotating neutron star, for arbitrary oscillations of a neutron star with arbi- 
trary magnetic field configuration under assumption of low current density in the inner 
parts of the magnetosphere. As an application we consider toroidal oscillation of a neu- 
tron star with dipole magnetic field and calculate energy losses. For some oscillation 
modes the longitudinal electric field can not be canceled by putting charged particles 
in the magnetosphere without a presence of strong electric current (j ~ -jj- p G3 c) . It 
is shown that the energy losses are strongly affected by plasma in the magnetosphere, 
and cannot be described by vacuum formulas. 

Key words: starsmeutron, oscillations, magnetic fields - pulsars: general 



INTRODUCTION 



Oscillations of neutron stars could give an opportunity to study their internal structure. Oscillations of a neutron star (NS) 
may be observable in pulsars just after a glitch. The energy release during the glitch is estimated as 



E ~ - 



2 M NS 



(1) 



where / is the moment of inertia of the NS, Q is an angular velocity of rotation, AQ, is a jump of angular velocity, Afj^-g is a 
mass of the NS, M CIUS ^ is a mass of the crust. For the Vela pulsar E ~ 10 37 ergs. There are about 20 glitching pulsars, which 
had suffered a total of 45 glitches (Lyne 1996). The energy release in each of these glitches is similar to that of the Vela pulsar. 
If several percent of this energy goes into excitation of oscillations it would be possible to observe them. Another example how 
oscillations of NS could be observed in pulsars is by observations of microstructure in single pulses. Boriakoff ( 1976] ) proposed 
that vibrations of a neutron star may cause the periodicity of micropulses observed in some pulsars. A difficulty in this theory 



is that no effective mechanism has been proposed for excitation of oscillations ( Hankins 1996 ) . 

Oscillating neutron stars have been proposed as a source for Galactic Gamma-Ray Bursts (GRBs) by Pachini & Rud- 



erman (1974) and Tsygan (1975). This idea is further developed by Blaes et al. (198!:), Smith & Epstein (1993), Fatuzzo 



& Melia (1993). Oscillation induced hard gamma radiation from NS was proposed by Bisnovatyi-Kogan (1995) and Cheng 



& Ding (1997) for explanation of the hard delayed emission observed from some GRBs (Hurley et al. 1994), even if the 



Gamma- Ray Burst itself is generated by a different mechanism (Bisnovatyi-Kogan et al. 1975). Recent GRBs observations 



make a Galactic origin of GRBs unlikely, but a Galactic model may not be excluded completely. At least for Soft Gamma 
Repeaters, which are believed to be neutron stars, oscillations of star can play an essential role (Duncan 1 998). 



1985 



Carroll et al. 1986 



Eigenfrequencies and eigenfunctions of neutron star oscillations have been computed by many authors ( McDermott et al 
[Btrohmayer 1991 ). These calculations have shown that typical periods of neutron star oscillations 



range from 0.1 up to several tens of milliseconds. On the other hand, observed pulsar radiation in different spectral regions is 
generated mostly in their magnetospheres (Triimper & Becker 1995). If we are looking for oscillation induced radiation from 
pulsars we should investigate processes in the magnetosphere. A complete solution for vacuum electromagnetic fields near an 
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oscillating magnetized neutron star was obtained by Muslimov and Tsygan ( 1986j ) . But for typical NS's magnetic field strength 
(> 10 9 G) the electric field arising from star's oscillation will be strong enough to pull charged particles from the surface in 
the magnetosphere. Hence, any realistic model of the magnetosphere of oscillating NS must take into account the presence 
of charged particles. This will also affect the electromagnetic energy losses of a pulsating magnetized star. Electromagnetic 
energy losses of an oscillating NS in vacuum were calculated by McDermott at al. (1984) and Muslimov and Tsygan (1986) 

As a first step in the investigation of magnetosphere of an oscillating neutron star with nonzero charge density we consider 
the inner part of the magnetosphere, assuming low current density there. We generalize a formula for Goldreich- Julian charge 
density, derived originally for rotating NS, for a NS oscillating in any arbitrary mode. This allow us to take into account 
plasma present in the magnetosphere of the NS and calculate the energy losses. The plan of the paper is following. In the 
first section we introduce some definitions and give algorithm for calculation of Goldreich- Julian charge density for arbitrary 
oscillation mode and magnetic field structure under the assumption of low current density, and discuss the limitations of this 
assumption. In the second section we apply this formalism to toroidal oscillations of a NS with a dipole magnetic field. We 
calculate the Goldreich- Julian charge density and energy losses for some oscillation modes and discuss factors affecting energy 
losses of oscillating neutron star. 



2 GOLDREICH- JULIAN CHARGE DENSITY FOR SMALL AMPLITUDE OSCILLATIONS OF A 
NEUTRON STAR. GENERAL FORMALISM. 

2.1 Basic definitions. 



As it was firstly pointed by Goldreich and Julian ( Goldreich fc Julian 196S ) a rotating NS can not be surrounded by vacuum. 
The electric field generated by rotation of magnetized NS will pull charged particles into the magnetosphere - the electrostatic 
force near the NS's surface in vacuu m would be much stronger than the gravitational one for both electrons and ions. According 
to the calculations of Jones (1986), the binding of charged particles in the crust of a typical pulsar is not strong enough to 
prevent them being pulled into the magnetosphere by the vacuum electric field. But even if the binding of charged particles in 
the cold crust of an old NS can prevent them escaping, the magnetosphere of a strong magnetized neutron star can be filled 
by charged particles by the mechanism proposed by Ruderman and Sutherland (1975) for pulsars. 

On the NS surface the vacuum electric field generated by rotation or oscillation has a radial component, whose value is a 
substantional part of the full field strength (Muslimov & Tsygan 1986; Goldreich & Julian 1969). To order of magnitude this 
field for the rotating star is 



E 1 - 



fir. 



B 



(2) 



where Q is the rotational frequency of the NS, r* is NS radius and B is the magnetic field strength near the star. In the case 
of NS oscillations the corresponding vacuum electric field is 



E° 



^B 



(3) 



where uj is the oscillation frequency, £ is the displacement amplitude. The vacuum electric field strength near the NS oscillating 
with period r will be of the same order as the field strength generated by rotation of the NS with period T e fr: 

T cB ~ r r j . (4) 

For typical oscillation parameters (McDermott et al. 1986) the strength of the vacuum electric field of the oscillating NS 
will be the same as that generated by a ~ 1 sec pulsar. Hence, the magnetosphere of an oscillating, even non-rotating, NS 
should be filled by charged particles. In the strong magnetic field of the NS [B > 10 9 G) synchrotron energy losses of charged 
particles in the magnetosphere are very high. They lose their perpendicular momentum (if any) very rapidly and occupy the 
first Landau level, i.e. they can move only along magnetic field lines. 

In the presence of strong longitudinal (parallel to B) electric field E\\ charged particles are accelerated to high energies 
and their curvature radiation in strong magnetic field produces an electron-positron pair cascade (Sturrock 1971; Daugherty 
fe Harding 1982). The particles produced in the cascade cancel the accelerating electric field. Because of this there should 



be a regular force-free (En <C B) configuration of the electromagnetic field, at least in the inner parts of the magnetosphere, 
where the NS magnetic field is strong enough to allow single photon pair creation: 



B = 0. 



(5) 



Let us call the electric field E G1 the Goldreich- Julian (GJ) electric field. We introduce the generalized Goldreich- Julian charge 
density as the charge density in the magnetosphere corresponding to the GJ electric field E : 



E" 



(6) 
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It follows from the above discussion, that the charge density in the inner parts of the magnetosphere of an oscillating NS 
should be approximately equal to the GJ charge density. 

In the following we are looking for the GJ electric field and the GJ charge density, because they should adequately 
describe the electric field and the charge density in the inner parts of the NS magnetosphere. Knowledge of them will allow 
us to calculate electromagnetic energy losses of an oscillating NS. 



2.2 Basic equations. Low current density approximation. 

We assume the neutron star to be a magnetized conducting sphere of a radius r„. We are interested only in the pulsation 
modes with non vanishing amplitude £ on the surface. Consider a region close to the NS (near zone), at distances from the 
NS surface smaller than the wave length A = 2irc/ui, where ui is the oscillation frequency and c is the speed of light. In the 
near zone r < A one can neglect the displacement current term. Maxwell's equations for the electric and magnetic fields in 
the near zone are: 

V E = 4tt P , (7) 
VxE = --d t B, (8) 

V B = 0, (9) 

4-7T 

VxB = —j, (10) 

c 

where B and E are the magnetic and electric fields, p and j are the charge and current density. We use in the notation 
for partial derivatives: d a = d/da. On the unperturbed surface of the NS E and B must satisfy the following boundary 
conditions: 

B r (r.) = B 0r , (11) 

47T 

Be,4,(r*) — Boe,o0 ± — J<t>,e , (12) 
EeA r *) = --(«xB )w, (13) 

E r {r,) = --(vxB ) r +4ira , (14) 
c 

where the subscripts r, 9, <j> denote vector components in a spherical coordinate system, v is the velocity of oscillation of the 
NS surface, Bo is the surface magnetic field inside the NS, a and J are the induced surface charge and current density. Here 
([H]) and (^|) are used to determinate the surface charge and current densities. Close to the NS the current flows along the 
magnetic field lines. So, for the inner parts of the magnetosphere it can be expressed as 

j = a(r,6,4>)B, (15) 

where a is a scalar function. This system must be completed by equation (||), defining the GJ electric field. We failed to find 
an analytical solution of the system in general case. But under some physical assumptions it can be solved analytically for an 
arbitrary oscillation mode and magnetic field configuration of the NS. 

Let us assume that the physical current density in the magnetosphere is low enough - the magnetic field to first order in 
£ = £/r* can be considered as generated only by volume currents inside the NS and by surface currents on its surface, i.e. 

-i«Vxfl (1) , (16) 
c 

where B^ is the first order in | perturbation of the magnetic field. To order of magnitude this means: 

i«i(^) c= ^i c (^), ( 17 ) 
r V r*/ cr» \LorJ 

where B^ is the unperturbed magnetic field strength. The value of the Goldreich- Julian charge density near the surface of 
the NS is, to order of magnitude, pG.i{r*) — B ^ V , where V is the oscillation velocity amplitude. Thus equation ( |l6| ) implies 

J<*H(r.)c (-£-) . (18) 

We call assumption ( |l6|) the low current density approximation. For regions of complete charge separation, where there are 
charged particles of only one sign, the maximum current density is pc.i c. The absolute value of the Goldreich- Julian charge 
density decreases with increasing r and in the near zone r <§C c/lu. Consequently condition ( |l^ ) for a charge-separated solution 
is satisfied automatically. Charged particles in the near zone move along magnetic field lines. Consequently the current flowing 
through a field line tube remains the same. Because of this, condition (|l8|) is automatically satisfied along field lines which 
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have crossed a region of charge separation. For regions in the near zone, where there are charged particles of different sign 
and the magnetic field lines have not crossed a charge separated region, condition ( |li| ) can be violated. 

We assume condition ( |l8| ) is satisfied in the whole near zone and find the GJ electric field and the GJ charge density. For 



some oscillation modes pcj obtained under this assumption has singularities. Because of the reasons discussed in Section 2.2 
a regular solution of the system (|| |7|jl5|) should exist for any oscillation mode and unperturbed magnetic field configuration 
of a NS. Hence, in cases where our solution has singularities, the low current density approximation fails and the physical 
current can not be neglected in the whole near zone. In some regions the current density will be to order of magnitude 

i^W'(^) • ( 19 ) 

This situation will be considered in subsequent papers. As we will show, the small current density approximation holds in 
regions of open field lines. In the whole near zone it is valid for more than 50% of the modes, at least for toroidal oscillation 
of a NS with a dipole magnetic field. 



2.3 Goldreich-Julian charge density. 

2.3.1 Equations for Goldreich-Julian charge density 

In low current density approximation one can neglect current term in Ampere's low (0). To the first order of dimensionless 
oscillation amplitude £ in the near zone we have: 

VxB = 0. 

Using the properties of solenoidal vector fields we write B in the form 
B = V X V X (Pe r ) + V X (Qe r ) , 

where e r = r/r. Scalar functions P and Q can be expanded in spherical harmonics Yi m as 



P(r,t) = 
Q(r,t) = 



E 

in i 

E 

im 



P m (r,t) Y lm (6,. 
Qim(r,t) Y lm (e, 



Substituting B in the form (|^) into the equation ( poj ) and multiplying the result by e r we get 
A n Q = 0, 

where An is the angular part of the Laplacian, 

1 1 
An = -— -dg(sm0 do) + . 9 . d&s . 
sin sur 

Substitution of the expansion (^3|) into the equation ( pi) ) gives us Qi m = for each l,m. Hence, Q = (see (Muslimov & 
Tsygan 1986)) and the magnetic field can be expressed in terms of only one scalar function P as 



(20) 
(21) 

(22) 
(23) 

(24) 
(25) 



B = VxVx(Pe r ) . 
Substituting B from expression 
1. 



fj in the Faraday's law Mj we get 



-VxVx(d t Pe r ) 



Integrating this equation: 
E = -ivx(<9 f Pe r ) - V* , 



(26) 
(27) 
(28) 



where "J is an arbitrary scalar function. 

From the theory of partial differential equations ( Elsgolts 1965| ; Brandt 1947) it is known that for any vector field A there 
exist a vector field perpendicular to A if and only if 



A ■ (VxA) 



(29) 



From equation ( pc| ) it follows that magnetic B satisfies equation (|2ij) and there always exist a vector field perpendicular to 
B. Let us assume that there exists an electric field E GJ with no component parallel to B satisfying the boundary conditions 
([t3|),(|l4|). Evidently it satisfies Maxwell equations (0)-(|2^) and consequently has the form 



E ( - 



= --Vx(d t Pe r ; 



(30) 
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The electric field E GJ satisfies equation (^). Substituting E GJ from the expression ( po[ ) in the equality (^) we have an equation 
for *gj 

-Vx(d t Pe r )-B + V* G j--B = 0. (31) 



c 



Substituting iS GJ from the expression (^o|) into equality (^|) we get an expression for the GJ charge density in terms of GJ 
potential ^gj'- 

Pgj = — — — A*&gj • (32) 

47T 

In other words, if one uses representation (^) for a magnetic field in the near zone, then the Goldreich- Julian electric field is 
written as a sum of two terms (^). The first one is the vacuum term, and the second one (— V^gj) represents the contribution 
of the charged particles in the magnetosphere of the NS. The potential is a solution of the equation ([jl|). In spherical 
coordinates, equation ( [3l| ) becomes (using expression |Pq| ): 

A n P cV* GJ - d r deP 9 e * GJ - -J^- d^P 9 * GJ + (d^P d e d t P - d r d e P d^dtP) = . (33) 

sir 9 c sin 9 

Let us consider small oscillations of a NS, ( < 1. We expand the function P in a series in £ and approximate it by the 
sum of the first two terms: 

P(t, r, 6, </>) « P (r, 6, <t>) + SP(t, r, 9, <t>) , (34) 

where the function Po{r, 9, <j>) is responsible for the unperturbed magnetic field and 8P(t, r, 9, <j)) is the first order term in an 
expansion of Pit, r, 9, <j>) in £. We expand SP in series of spherical harmonics 

6P = Y,SPi™(t,r)Y lm (e,<P). (35) 

Subsequently we will need only the time-derivative of SP. After substitution of the expansion (^) in the equation ( |2c| ) we 
find (see Appendix [X|) 

dt SP = (^) 1 d t Sp lm (t) Y lm (9, 4>) ■ (36) 

l,m 

The coefficients dtSpi m (t) are fixed by the boundary conditions. 
2.3.2 Boundary conditions 

The boundary conditions for the potential "Jgj are obtained from the boundary conditions for electric field (|l3|). The tangential 
components of the GJ electric field on the surface of the NS are 



Ef 3 \ r =r, = ~~ (-AnP Va, + drdd,P Vr 

c r \r sm 9 



— - (-AnP Vg + d T dgP V, 

c r \r 



(37) 
(38) 

From the expression ( pfj| ) the tangential components of the GJ electric field outside the NS are 

Ef " = - i fi J- d^dtP + de* G j) (39) 

= - (-dedtP - -^r d^J) ■ (40) 
r \c sm9 J 



Equating these expressions on the surface of the NS and eliminating v&gj one obtains boundary conditions for the 9— and . 
derivatives of $ G .i 



9e , I / Gj|r=r», = 



-- ( -AnP vs + -r—; d r dsP v r + d<hdtP 
c \r smfc* smfl 

9</, , l V Gj|r=r, = — (- sin9 AnP vg + sm9d r deP v r + sm9dedtP 



(41) 



(42) 



From these boundary conditions we obtain the boundary condition to be applied to Spi m (t). Differentiating equation ( pfl| ) 
with respect to </> and equation (^) with respect to 9 and equating results we get an expression for 5pi m (t) to first order in f 
(see Appendix ^) 

d t S Plm (t) = ^ * ^ J dnY t * m [v V(AnPo)+A n P (V ■v ± )+r 2 (V±(d r P )- V x )«r]| r=r< . (43) 
© 0000 RAS, MNRAS 000, 000-000 
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Boundary condition for ^g.i can be obtained by integrating equation ( f4l| ) or equation (^). For convenience we will use as 
boundary condition the result of integrating equation ( fi"l| ) over 6. For a perturbation depending on time t as e -1 "* we get to 
first order in £ 



*G.,|r=r» =-- d6 (-AnP v<j, + -— drd^Po v r + -r-z d<pd t 5P 
c J \r smO sine* 



e-' Mt F(<t>) , (44) 



where F is a function of 4>. 

For each oscillation mode the corresponding velocity field is continuously differentiable. From the boundary condition for 
the electric fiel d (|13| ) it follows that the tangential components of E GJ are finite. The vacuum term on the the right hand side 
in expression (BOl) for the electric field E GZ is also finite (with the natural assumption that P is continuously differentiable 
with respect to 8 and (f> expressions (^3|) and ( pj| ) are finite for r = r»). Let us consider the azimuthal component of GJ electric 
field E GJ near the polesQ In expression ( [io| ) for E GJ outside the NS both E GJ and the first (vacuum term) on the right side 

is also finite, hence d^G]\{e=o,n; r=r»} = 0. Thus 



of expression (^), are finite. Consequently the second term ( — 

r— r* 

in the boundary condition ( p4| ) one must choose F such that Vl?Gj|fl=o,7r; r=r» 
Using the gauge freedom we choose 

*G.l|fl=0;r = r» = . (45) 

Using approximation ( [34] ) for P we write an equation for \&gj in the case of small oscillations 

AnP a r * GJ - cWePo c*e* G j - — !~-;<9r^Po 9*¥gj + — ^ d B d t 8P - d r d e P d^dtSP) = . (46) 

sur 6< csmv 

The first-order partial differential equation ( firjl ) together with the boundary conditions (^), (^) and dt&P described by 
formulas (^6|), ( fUjl ) determines \I/gj. The solution of this problem provides ^gj to first order in f. The Goldreich- Julian charge 
density p G ,i can then be obtained from the formula 



2.4 Rotating NS - Pulsar 



At the end of this section we consider an important particular case - a rotating neutron star. We choose 2-axis to be parallel 
to the rotation axis. In this case the partial derivative dt can be replaced by — fldj,, where Q is the angular velocity. Thus 



Equation ( |33| ) for ^ g .t takes the form 

1 



A n P ft.* GJ - drdeP %* GJ 
By direct substitution of 



0*GJ - 



n 1 



c sin 1 



{drd^P ded^P - d r deP d^P) = . 



■ sin 9 d e P 



(47) 
(48) 
(49) 



it can be shown that the potential (^) is a solution of the equation (fig), and satisfies the boundary conditions ( fi"l| , ^). We 
note, that equation ( |33| ) and its particular form (j48|), and boundary conditions (^l|, ^) are valid for an arbitrary oscillation 
amplitude, i.e. also for rotation of the NS. Substituting potential ^/^j into the expression for the GJ electric field ( |30| ) one 
gets 



-(nxr)xB . 



(50) 



This result was obtained by Goldreich and Julian (1969), and Mestel (1971). 



3 GOLDREICH- JULIAN CHARGE DENSITY AND ELECTROMAGNETIC ENERGY LOSSES 
OF A NEUTRON STAR WITH DIPOLE MAGNETIC FIELD. SMALL AMPLITUDE TOROIDAL 
OSCILLATION. 



points with 8 = 0(it), 
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3.1 General Formulas 

As an application of the developed formalism we consider the case of small amplitude toroidal oscillations of a NS with dipole 
magnetic field. Velocity field on the NS's surface for a toroidal oscillation mode (l,m) is described by (Unno et al. 197£ ) : 



0, Vg = e 



W(r) — d^Y lm {e,4>), v+ 



- iu3t W(r)dgY lm (9^) 



(51) 



where W is transversal velocity amplitude. For simplicity we assume that the mode axis is aligned with the dipole moment fi. 
We can do this without loss of generality because any oscillation mode with mode axis not aligned with the dipole moment 
can be represented by series of oscillation modes with mode axis parallel to fi. In this case the unperturbed magnetic field is 



B = B, 



■ (r±y cos e er + ±B (^y 



sm U eg 



(52) 



where e r and ee are unit coordinate vectors. This field is described by the scalar function Pq' p (see eq. [Q, J34|) according 
to the formulas (A7), (A£) : 



(53) 



Scalar function <5P dlp (see eq. |2q]) describes first order inf magnetic field perturbation. The time derivative of this function 
for an oscillation mode (l,m), according formulas (^), (B7), is 



dtSP^ P = B i 



+ 1) ( r ) 



1(1 + 1) 



(54) 



where w = e l0Jt W. Substituting P^ lp and dtSPj^ into the general equation ( pi) ) we get a partial differential equation for GJ 
potential ^Jjf for a NS with dipole magnetic field oscillating with a small amplitude in a toroidal mode (l,m): 

m 2 B w /r,\ 1+1 



2 cos 9 drV l £ + -sm6 d fl tf - — ■ 

r 1(1 + 1) c 



Yh 



0. 



(55) 



where corresponds to one excited oscillation mode (l,m), and for general mode ^gj = Xw m *gj- The characteristics of 
the equation (m) are 



(t) 



Bowr, fr, 



1(1 - 



Y lri 



t 


= Co 




= Ci 


1/2 






= c 2 


dO 


= c 3 



(56) 



(sin60 2i+1 

The integral of the equation ^Ej) is an arbitrary function of constants Co, Ci, C2, C3 

ip(C ,C l ,C a ,Ca) = Q, 



Expressing we have for the general solution of equation (|55| 



Bowr* ( r* 



1(1 + 1) c 



(7)' 



Y lri 



(sin t 



■ de + 



sin 9 



1/2 



(57) 
(58) 



where <E> is an arbitrary function. In order for represented by the expression ( p&j ) to be a GJ potential it must satisfy 
the boundary conditions on the surface of the NS (Q) and (0), which in the case of toroidal oscillations of a NS with dipole 
magnetic field take the form 



\Si lm \ 
*6J r=r, 



^GJ \e=0,r=r, 



Bowr f 



COS 6 dfiYl, 



Y, r 



1(1 + 1) sin 6 



d9 + e~ iu,t F(<j>) , 



. 



(59) 
(60) 



Substituting given by the expression (g8|) into the boundary conditions (g9|), ( |60| ) we get the boundary condition for the 
function $ im : 

™2 



$ im I sing 



1/2 



Bowr t 



cosOdgYh 



Y lr 



Bowr* m . 2 i , 
■ sm ( 



c 1(1 + 1) 
Bowr t 



1(1 + 1) sin 

Ylm 



de 



cos 9 dgYi r 



(sin(9) 2 ' +1 
1(1 + 1) sir 



d9 + 



AO 



(61) 
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where the last term is added in order to satisfy the second boundary condition (|60]). The function & lm depends on 9 in the 
combination ' (77) 12 sin To get & lm one has to express the right hand side of equation (j6lj) in terms of 'sin #' and to replace 

'sin#' by ' (^r) ^ 2 sin 9'. Substituting the function «J >!m into the expression we get the potential "^^j for small toroidal 
oscillations of a NS with a dipole magnetic field for the mode (2, m). Using the formula (^2|) we get the Goldreich- Julian 
charge density for that mode. We have developed a set of programs on the computer algebra language MATHEMATICA 3.0 
(Wolfram 1996) for calculating the analytical expressions of Pg"j according to the algorithm described in this section. These 
programs were tested for some (I, m) by comparing results obtained by hand and by computer, and also by checking the 
condition E GJ ■ B = 0. 



3.2 Main Results 

3.2.1 Goldreich-, Julian charge density 

In the solution for the whole region the expression ( |6l| ) is used, where on the right hand side 'cos 8' should be replaced by 
'(l — sin 2 0) ^ ' f° r < 9 < t; (first hemisphere) and by '— (l — (^) sin 2 ^ ' for ^ < 9 < ty (second hemisphere). 
Because of this there are two different expressions for ^gj for both hemispheres. If these expressions give different results 
for 0=2 for r > r* there is a discontinuity in the function ^gj at the equatorial plane, and consequently pci(r, ? ,cj>) 
becomes infinite. As we discussed in section this unphysical result indicate that low current density approximation for 
such oscillation modes cannot be applied in the whole near zone. In this situation it is impossible to cancel accelerating electric 
field in the whole near zone without presence of a strong electric current (|l^) in some magnetospheric regions. Examples of 
such oscillation modes are modes (I, m) = (2, 0), (3, 3). The corresponding velocity fields are given in Fig.[j]. The dependence of 
the potential 'J/ 1 ^ for these modes along dipole magnetic field lines is shown in Fig^||| In the following we normalize the G-J 
potential *P G .i and charge density p G j by BoWr,/c and BoW/r*c respectively. The jump of $ GJ (r, 7r/2, (/>) decreases, as the 
angle at which corresponding magnetic field line intersect the NS's surface is increasing, ^gj is a continuous function of 9 on 
the surface of the NS, nevertheless GJ charge density diverges on the equator even on the surface of NS (see Fig|B|, In these 
figures the charge density pGj|r=r, is shown in spherical coordinates ( PGj|r=r, | ,9,cj>). Here the radial coordinate represents 
absolute value of the GJ charge density as function of the polar angle 9 and the azimuthal angle cj>. Analytical expressions for 
the GJ charge density for the discussed modes are given in Appendix ^. At the equatorial plane p^i i s infinite. The mode 
(2, 1) is a representative of another class of oscillation modes, in which ^(ji is a continuous function of 9 (see Fig.^) and p^ 
is finite everywhere. The velocity field and p G j are shown in FigsjjJ |^ For such modes it is possible to cancel longitudinal 
electric field without generating strong current along magnetic field lines. 

Consider now the difference between two classes of toroidal oscillation modes for a NS with a dipole magnetic field. The 
boundary condition for is proportional to the vector product of v and B. The axis z is directed along the dipole magnetic 
momentum. If the velocity field is symmetric relative to the equatorial plane (I — m is an odd number), then the boundary 
condition for i s a l so symmetric relatively to the equatorial plane. Hence, after expression of all trigonometric functions in 
( |si| ) through 'cos 9' and 'sin 9' , 'cos 9 l appears in the right hand side of ( |6l| ) only in even powers. Consequently the function 
$ !m will be the same in both hemispheres and pjjj w iU De finite. If velocity field is antisymmetric relative to the equatorial 
plane (I — m is an even number or zero), then 'cos 9 l is contained in the right side of (|6TJ) in odd powers, and the function <l?' m 
may have discontinuity in the plane 9 = |-, and Pg"j f° r such modes becomes infinite at the equatorial plane. An example of 
oscillation mode with velocity field antisymmetric relative to the equatorial plane and for which ^[jr is continuous function 
of 9 is the mode (1, 1). Graphics for the potential *1/gj! Goldreich- Julian charge density p G j are shown in Figs. [| [| So, for a 
dipole magnetic field the longitudinal electric field generated by the oscillation modes with odd (I — m) can be canceled by 
putting of charged particles in the magnetosphere. Among the modes with even or zero (I — m) there are modes, for which 
longitudinal electric field cannot be canceled by putting of charged particles without presence of strong electric current along 
some magnetic field lines. 

For an oscillation mode with odd (I — m) there is a regular solution for p G j and E GJ in the low current density approx- 
imation. This solution will be stable because of Lenz's low: increasing of the current will lead to a generation of a magnetic 
field inducing an electric field, which will prevent the current to grow. In other words, configuration with the smallest possible 
current will take place. For oscillation modes with even or zero (I — m), where no regular solution for GJ charge density and 
electric field exists, the low current density approximation cannot be used in the whole near zone. Because some of oscillation 
modes of this class (at least one) possess regular solution in low current density approximation, the total amount of modes, 
for which this approximation may be used in the whole near zone exceeds 50%. 

The Goldreich- Julian charge density gives only a characteristic charge density in the magnetosphere. The particle density 
can be obtained by so lving full sy stem of MHD equations, but this has not been done even for an aligned rotator with a 



dipole magnetic field ( Michel 1991 ). The case of nonradial oscillation is much more complicated, because of the absence of 
stationarity and axial symmetry. About particle number density in the region of closed magnetic field lines we can say, that 
for oscillation modes with odd (I — m) there can be charge separated regions in the near zone. For such modes the foot points 
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of magnetic field lines have the same sign of the potential and there are regions where GJ charge density does not change the 
sign along the whole magnetic filed line. Hence, charged particles of only one sign can be pulled from the NS surface into these 
regions. Evidently for modes with even or zero (I — m) there can not be charge-separated regions in the near zone, because 
GJ charge density change the sign along any magnetic field line. 



3.2.2 Energy Losses 

In previous works the calculation of electromagnetic energy losses of an oscillating neutron star includes only radiation of 
electromagnetic waves in the vacuum. If there is plasma in the magnetosphere, then the energy will be lost by a transformation 
of the oscillation energy into kinetic energy of an outflowing plasma, as it was proposed for pulsar energy losses by Goldreich 



and Julian (1969). They obtained the energy losses of a rotating aligned NS through outflow of the charged particles from 
a region of open field lines, and found it to be equal to the loss in vacuum through radiation of electromagnetic waves by 
a perpendicular rotator. In this paper we consider electric and magnetic fields in the near zone, hence we can not explicitly 
show the existence of plasma outflow as in the aligned rotator by Goldreich and Julian, but a qualitatively the picture is the 
following. Electric current arising from the motion of the charged particles in the magnetosphere generate a magnetic field. In 
the region where the value of this field becomes larger than the unperturbed NS magnetic field the field lines become open, on 
the Alfvenic surface. Plasma escaping from the region of open field lines produces an electromagnetically driven stellar wind. 
The wind causes a net electric current which closes at infinity. This current flows beneath the stellar surface between positive 
and negative emission regions. Because it must cross magnetic field lines there, it exerts a braking torque on the oscillating 



NS and reduces the oscillational energy. This picture is similar to the one proposed by Ruderman & Sutherland ( 1975) for 
pulsars, but now we must get the boundary of the open field line region self-consistently. In our case we can not approximate 
the closed field line region by a "zone of corotation", rather we should self-consistently determine the last closed field line as 
the last field line lying inside the Alfven surface. 

Consider region near the pole. Let us denote as do a polar angle by which the last closed field line intersects the surface 
of the NS. Similar to pulsars above a polar cap an acceleration zone and a zone of e + — e~ pair generation above it will be 
builded. Electron-positron pairs cancel accelerating electric field, i.e. motion of charged particle above the accelerating zone 
will be not influenced by electric filed generated by stellar oscillations. In the region of open field lines, where particles escape 
from the neutron star and form relativistic wind, the average time of crossing of the acceleration zone by the particle is much 
less then the oscillation period (see Discussion), so particles practically do not return back to the star. Hence, averaged over 
the oscillation period, the energy loss through the outflow of plasma from open field line region is 



dt / \j {r^6,4>) A lm (e,cj))\ rism8A6A4>, (62) 

Jo Jo 

where A(9, <f>) is the work done by the electric field to move a unit charge to the point with coordinate (r„, 9, (ft): 

A(0,<f,)~ f Ef\r„8'^)r,A6' . (63) 



Because of reasons discussed at the end of section 2.2 charge density in the inner parts of the magnetosphere must be 
approximately equal to the Goldreich- Julian charge density p GJ . In the region of plasma outflow charged particles are streaming 
along magnetic field lines in the same direction with ultrarelativistic velocities. Charge density of outflowing plasma is equal to 
the Goldreich- Julian charge density. Hence, outflowing ultrarelativistic electron-positron plasma builds a net current density 

j k Pgj c . (64) 

The current density can differ from the values given by expression ( |64| ) due to difference in averaged velocities of electron 
and positron components of the plasma (discussion on this topic see in Lyubarskii ( 1992|)), but on the order of magnitude 



expression (|6j) gives good estimation for the current density (Ruderman & Sutherland 1975). Hence, for the open field line 
region condition ( |l8| ) is satisfied and in expression ( p^ ) we can use p GJ obtained by solving of equation (^(j) for any oscillation 
mode. 

The last closed field line is the line for which the kinetic energy density of outflowing plasma on equator (at the point 
(i? a , 7r/2, (j>)) becomes equal to the corresponding energy density of the NS magnetic field: 



AtyRIc 8tt 



(65) 

r T,r=Ra. 



Expressing the right hand side of this equation in terms of angle 8o, we get two equations (|62j), ( |65| ) for a self-consistent 
determination of do and . 
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We have solved these equations for oscillation modes (1,1), (2,0), and (3,0). Oscillation periods^] were taken from McDer- 
mott et al. (1988). The ratio of energy losses to the vacuum ones as a function of the amplitude of dimensionless transversal 
displacement £ for modes (1,1), (2,0), and (3,0) are given in Table El One can see that energy losses through outflowing 
plasma are much less than in vacuum for small values of £. On the other hand, for oscillation modes I > 2, for displacement 
amplitudes £ larger than a value |cr, energy losses exceed the vacuum losses (see Fig. [ilj). For some oscillation modes for 
which 6>o <C 1, we can linearize equations (^), ( pB] ) in do and get an analytical estimation of 6*0. For toroidal oscillation mode 
(I, m) the velocity amplitude V near the pole is of the order 



V ~ w 



9 m -\ m/0 
9, m = 



The electric field is of the order E GJ ~ BV/c. From formula (KS3) we have 



BWr t 



9 m , m/0 
9 2 , m = 



(66) 



(67) 



The Goldreich- Julian charge density is 

BW m 
Pgj ~ 

Substituting expressions (^) and (^) into equations ( p3| ) and (^) we get the energy loss through the plasma outflow 

along open field lines: 



(68) 



R2 2 fW\ 

e p i ~ S r»c I — I 



m / 
m = 



(69) 



The angle 8o is determined from the equation (|65|). After substitution of e p j into equation (|35|) we get 
m / 



W\5 



(70) 



m = 



We see that the angle So is small only for modes with m < 3. This is because of the small value of V in the polar region for 
modes with m > 3. For such modes a linear analysis is not possible and the case of large 6q needs additional investigation, 
which will be given elsewhere. The energy loss through plasma outflow is 



D 2 2 

e p i ~ B r» c 



m = 1,2 



m = 



(71) 



The vacuum energy loss for the oscillation mode (l,m), according to McDermott et al. (1984) is 
(^) 8 , l = l, m = Q 



£vac ~ B^rlc i 2 



21 + 2 



(72) 



, I = l,ra = 1; I > 2,m = 0.. .1 
From these formulas we get the oscillation amplitude £ for which the plasma energy loss is larger than the vacuum loss 

/ = 1,771 = 



£ >Cc 



(¥)' 



\ c 



1/3-1 



/ > 2,771 = 0; I > 1,771 = 1 



/ > 2,m = 2 



(73) 



From formula @ it follows that for modes (I = l,m = 0), (Z > 2,m = 0, 1) and (I > 4,m = 2) nonvacuum energy loss 
exceeds the vacuum loss even for small displacement amplitudes £, because for these modes £cr <IC 1. 



t According to McDermott et al. (1988) the dependence of the eigenfrequency on I is very weak, at least for small I, so we assumed the 
same oscillation periods for modes (2,0) and (3,0) 
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Mode 


^pl/^vac 




(1.1) 


0.02 £ 2 


0.9/f 


(2,0) 


77i 2 r 2 


0-2/5 


(3,0) 


9 ■ 10 5 £ 2 t 4 


0.09/1 



Table 1. Leading terms in expansion over £ of ratios of electromagnetic energy loss through plasma outflow to the vacuum one e p i/e V ac 
and distance in the equatorial plane of the last closed field line to the wave length Rr™ /A for toroidal oscillation modes (1, 1), (2, 0) and 
(3,0). r is period of oscillations in milliseconds. 



Electromagnetic waves radiated by the oscillating NS are screened by the plasma. For small £ the nonvacuum energy 
loss according to (^) is smaller than the vacuum loss because at the GJ density of particles the accelerated plasma cannot 
escape from the magnetosphere. The ratio of the equatorial radius for the last closed field line, to the wave length A are 
given in Table [j]. For mode (3, 0) this dependence is shown in Fig. [ill. We see that the last closed field line for small values of 
displacement amplitude lies deep inside the wave zone. In the region limited by the last closed field lines and a sphere of radius 
A (formal boundary of the wave zone) in this situation confined plasma waves are excited such as Alfven or magnetosonic wave 
by the oscillation of the NS, see Fig. |l2|. This will lead to accumulation of energy in this region in the form of plasma wave 
energy. When the density of this energy becomes larger than the magnetic energy, field lines becomes open. These processes 
will lead to a decrease of the equatorial radius of the last closed field lines, and hence, a increase of the energy loss. 

A second process which could lead to larger energy losses than obtained by the solution of the equation ( jj^ ) , is rotation 
of the neutron star. In this case the equatorial radius of the last closed field line is the radius of the light cylinder Rl = fi/c, 
and surface area for plasma outflow may be considerably increased. It is not clear, however, how the interaction between 
rotation an oscillation will influence the total energy losses. 



4 DISCUSSION 



Under the assumption of low current density we have developed a formalism describing the inner parts of the magnetosphere 
of an oscillating neutron star in the same way as in the case of rotation. This assumption is valid in the region of open field 
lines for any oscillation mode. For some oscillation modes the assumption about low current density everywhere in the inner 
magnetosphere leads to unphysical results - the Goldreich- Julian charge density becomes infinite in some points. For such 
oscillation modes the low current density approximation cannot be used in the whole region of closed magnetic field lines 
near the NS. The current density for these modes along some closed magnetic field lines will be of the order of magnitude 
j ~ — Pgj c. For a rotating neutron star our formalism yields the same values for the Goldreich-Julian charge density and 
electric field as obtained in works of Goldreich and Julian (1969) and Mestel (1971). 

We applied the general formalism to the case of toroidal oscillations of a neutron star with dipole magnetic field. We 
calculated the GJ charge density and showed, that for the case considered the assumption about low current density in 
the whole inner magnetosphere is valid for more than half of all modes. We calculated the electromagnetic energy losses 
of a NS with dipole magnetic field for some toroidal oscillation modes and showed, that the energy losses of an oscillating 
neutron star are strongly affected by the plasma present in its magnetosphere. Electromagnetic energy losses of an oscillating 
NS due to plasma outflow from the magnetosphere have in general another dependence on oscillation frequency u> and the 
dimensionless displacement amplitude £ than the energy losses due to radiation of electromagnetic waves in vacuum. This 
affects all previous calculations of the electromagnetic damping rate of NS oscillation. Our calculations give a lower limit on 



electromagnetic energy losses of an oscillating NS because of the reasons discussed at the end of section 3.2.2. The energy 
of a star oscillation depends on the dimensionless displacement amplitude as £ 2 . Energy losses through plasma outflow in 
contrast to energy losses through radiation of electromagnetic waves in vacuum have another dependence on the displacement 
amplitude. Consequently the electromagnetic energy losses of an oscillating NS in general depends on the oscillation amplitude 
£. To estimate the oscillation damping time for the considered oscillation modes one can use values given in McDermott et 
al. (1988), Table 6, multiplying by functions from Table |l| of this paper for a given displacement amplitude. Here we have 
restricted ourselves to a simple case of toroidal oscillations. Starquakes should excite a whole spectrum of oscillations, including 
p- and g- modes. The magnetosphere structure produced by these modes will be considered in a future paper. 

With the proposed formalism it is possible to apply theoretical models developed for pulsars to oscillat ing n eutron 
stars. For example, one ca n inv estigate the acceleration mechanism proposed by Scharlemann, Arons & Fawley ( 1978 ) or by 
Ruderman & Sutherland (1975) for oscillating neutron star. The inertial frame dragging mechanism proposed by Muslimov 



& Tsygan ( 1992 ) does not work here, because we consider a nonrotating star. Characteristic time of pair cascade formation 
is estimated as /ippf/c, where /ippf is the hight of pair formation front above the surface. In "Ruderman-Sutherland" model 
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minimum value of /ippf is the same as for a NS rotating with the angular frequency Q e g = w£, and for usually assumed 
parameters of NS oscillation is of order 10 4 — 10 5 cm. In "Arons-Scharlemann" like model the hight of pair formation front can 
be smaller then in a pulsar rotating with the angular frequency Q e g (several km), because in general for nonradial oscillation 
ratio pcj/B along magnetic field lines increase faster than in magnetosphere of a rotating NS, what leads to larger accelerating 
electric fields. Because the characteristic time of pair cascade formation is much shorter than the period of oscillation, one can 
consider the polar cap acceleration zone in any given moment of time as stationary. Problem of return current region for the 
case of oscillating NS as in the case of pulsar remains open. Similar to pulsars return current can flow along the last closed 
magnetic field line. 

The authors wish to thank A. I. Tsygan for helpful discussions. This work has been supported by Russian Foundation of 
Basic Research (RFFI) under grants 96-02-16553 and 99-02-18180. 
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Figure 1. Velocity field on a sphere for toroidal modes (1,1), (2,0), (2,1) and (3,3) are shown at the time t = 2irn/u>, where n is an 
integer, as projection on the meridional plane <f> = —115°. Circles corresponding to the longitudes 0° (right), -90° (left) and latitude 0° 
are shown by dashed lines. 
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8 

Figure 2. The potential ipQ, along a dipolar magnetic field line as a function of the polar angle 9 is shown for 5 field lines with azimuthal 
angle cf> = n/2, for t = 2-nn/u), at which the maximal absolute values of the potential are reached. The values of 9 of the left and right 
ends of the lines 1-5 determine the polar angle at which corresponding magnetic field line crosses the surface of the neutron star. 




Figure 3. left: The charge density PGj|r=r, ' s snown in a spherical coordinate system (|pc^j|r=r*| ,9,<j>) for t = 2-nn/u. The radial 
coordinate represents the absolute value of the GJ charge density as function of the polar angle 9 and the azimuthal angle cf>. Positive 
values of charge density are shown by a gray surface and negative ones by white. The dipole magnetic moment fi is directed upwards 
along the z axis. Azimuthal angle is counted from the x-axis, as in Fig. hi. right: Cross-section of surfaces from the left figure by the plane 
<f> = it/2 , where the maximal absolute values of the charge density P G ^ J |r=r» are reached. Circles correspond to the values 0.01, 0.02, 
and 0.03. Positive values are shown solid, negative ones dashed. 
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Figure 4. The potential ^f^Pj along a dipolar magnetic field line as a function of the polar angle 9 is shown for 5 field lines for t = 2irn/uj, 
at which the maximal absolute values of the potential are reached. The values of 9 of the left and right ends of the lines 1-5 determine 
the polar angle at which corresponding magnetic field line crosses the surface of the neutron star. The discontinuity is shown by the 
jump of if) corresponding to the same field line in the opposite hemispheres. 



Z 




Figure 5. left: The charge density p^?j\r=r„ is shown in a spherical coordinate system (|pc^}|r=r„| ,9,cf>) for t = 2-wn/u). The radial 
coordinate represents the absolute value of the GJ charge density as function of the polar angle 9 and the azimuthal angle cf>. Positive 
values of charge density are shown by a gray surface and negative ones by white. The dipole magnetic moment p is directed upwards 
along the z axis. Azimuthal angle is counted from the x-axis, as in Fig. |[| Note that for 8 = tt/2, p^P 3 is infinite, right: Cross-section of 
surfaces from the left figure by a meridional plane. Circles correspond to the values 0.1, 0.2, and 0.3. Positive values are shown solid, 
negative ones dashed. Note that for 9 = n/2, p^Pj is infinite. 
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Figure 6. The potential along a dipolar magnetic field line as a function of the polar angle 8 is shown for 5 field lines with azimuthal 
angle cf> = n/2, for t = 2irn/u), at which the maximal absolute values of the potential are reached. The values of 8 of the left and right 
ends of the lines 1-5 determine the polar angle at which corresponding magnetic field line crosses the surface of the neutron star. 





0.05 



Figure 7. left: The charge density PGj|r=r« is shown in a spherical coordinate system (|/»oj]r=)-. | ,8,4>) f° r * = 2-wn/u). The radial 
coordinate represents the absolute value of the GJ charge density as function of the polar angle 8 and the azimuthal angle cf>. Positive 
values of charge density are shown by a gray surface and negative ones by white. The dipole magnetic moment fi is directed upwards 
along the z axis. Azimuthal angle is counted from the x-axis, as in Fig. [l]. right Cross-section of surfaces from the left figure by the plane 
(j> = n/2 , where the maximal absolute values of the charge density reached. Circles correspond to the values 0.05, 0.1, and 

0.15. Positive values are shown solid, negative ones dashed. 
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Figure 8. The potential \Pqj along a dipolar magnetic field line as a function of the polar angle 8 is shown for 5 field lines with azimuthal 
angle cf> = 7r/6, for t = 2nn/u>, at which the maximal absolute values of the potential are reached. The values of 8 of the left and right 
ends of the lines 1-5 determine the polar angle at which corresponding magnetic field line crosses the surface of the neutron star. The 
discontinuity is shown by the jump of ip corresponding to the same field line in the opposite hemispheres. 




Figure 9. left: The charge density Pcfj!r=r» is shown in a spherical coordinate system (|p[^J| r — rt | ,8,<j>) for t = 2nn/u>. The radial 
coordinate represents the absolute value of the GJ charge density as function of the polar angle 8 and the azimuthal angle cf>. Positive 
values of charge density are shown by a gray surface and negative ones by white. The dipole magnetic moment fi is directed upwards 
along the z axis. Azimuthal angle is counted from the x-axis, as in Fig. hi. Note that for 9 = tt/2, p^ is infinite, right: Cross-section 
of surfaces from the left figure by the plane <f> = it/ 6 , where the maximal absolute values of the charge density p^j\ r =r, are reached. 
Circles correspond to the values 0.04, 0.08, and 0.12. Positive values are shown solid, negative ones dashed. Note that for 8 = n/2, p^?j 
is infinite. 
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Figure 10. Ratio of the energy loss through outflow of plasma from the region of open field lines to the loss in vacuum is shown as a 
function of the dimensionless displacement amplitude £ for the mode (3,0). The solid line corresponds the oscillation period T = 19ms 
(mode 3to in notations of McDermott et al. (1988)), the dashed line - T = 1.7ms (mode 3*1), the dotted line - T = 1.08ms (mode 3*2); 
see note on the page [To} 
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Figure 11. Ratio of the equatorial radius of the last closed field line to the wavelength A is shown as a function of the dimensionless 
displacement amplitude £ for the mode (3,0). This ratio does not depend on the oscillation frequency w. 
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Figure 12. Structure of the magnetosphere of an oscillating neutron star (arbitrary scaled). See explanations in text. 
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APPENDIX 
APPENDIX A: CALCULATION OF P AND SP 

Substituting B from formula ( |26| ) using expansion into the equation ( |2o| ) we have 
V X V X V X (P + 8P)e r = . 



Because P does not depend on t, equation (Al) is equivalent to the following two equations: 

Vx Vx Vx(P e r ) = 
VxVxVx(<5Pe r ) = . 



In spherical coordinate equation (A3) is 



VxVxVx(<5Pe r ) =- ee —L- (J-d <l) A n SP + d rr d < t,Sp)+e <l) - ( — d e A n SP + d rr d e SP) = 0. 

rsmO \r 2 J r \r 2 J 



Substituting expansion of SP (|3jj) we get 

+ l)Splm(r,t) + drr8pl m (r,t) = . 



The solution of equation flA5|) which vanishes at infinity is 
8pi m (r,t) = (j^) Spi m (t) . 

The time-dependent coefficients p; m (i) have to be determined from the boundary conditions. Substituting 
into expansion ( ^ ) and taking the derivative with respect to t we get expression (^6|). 
Similarly, for the coefficients Poim in the expansion of the function Po we have 

Poi m (r) = pim 

Substituting the expansion of the function Po in spherical harmonics into the boundary condition (pl|), 
expression by Y* m , and integrating it over solid angle 4ir, we have (compare with formula (16) in Muslimov < 



(Al) 

(A2) 
(A3) 

(A4) 
(A5) 
(A6) 

expression (A6) 
(A7) 

multiplying this 
i Tsygan ( |l986| )) 



Plm = 



rl 



1(1 + 1) 



Bo r Yi m dfi, 



(A8) 



where Yi m is inverse spherical harmonic J Yi m Y^ m , dfi — 6n>6 mm i and dfi = sin9d9d(f> is the solid angle 



APPENDIX B: CALCULATION OF 6p lm (t) 

Differentiating equation ( |il| ) with respect to <j> and equation (^) with respect to 9, equating the results and multiplying by 
1 / sin 6 we have 



±-d B {sm6 de{d t P)) + -L-d^dtP) 
in 9 sm z 6 



A n P 
1 



-J—. -d e (sm9 v e ) + -^-^ds v& 
smO suit) 



d e (sin OdrdeP) + 



sm 

-deAnP vg H —. - 

r r sm 9 



-drdj,j>P 



+ 

V r + 



d<p AnP v<p 



+ 



d r dgP dgv r + 



-drdcpP d<pV r 



(Bl) 



Simplifying equation (Bl) by writing P as a sum of the two terms (|34j), we get an expression for time-dependent part SP to 
first order in £ 



d t Aci6P\, 



(v r drAnPo + vg -dgAnPo + v^ —^-rd^AnPo) + A n P ( —^-:dg(sm 9 v g ) + — 
\ r r sm 6 J \ r sin 9 rsm9 

(drdgPo d e v r + . \ d r d^P d^v r \ 
\ sm 9 J 



or in vector notation 

d t A Q 8P\ r=rr = - {v ■ V(AnPo) + A«P (V ■ v±) + r 2 (V ± (d r P ) • V_l) w r }| 



(B2) 
(B3) 
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where v± is the tangential part of the velocity, v± = eg vg + and Vi is the angular part of the gradient: Vi = 

eg l/r dg + e^ l/(rsin#) 8$. Substituting the expansion of SP in spherical harmonics (pr]]), multiplying the result by Y{ m 
and integrating it over the solid angle 4ty, we get an expression for the coefficients of dtSp(t) (^) in the expansion of SP in 
spherical harmonics: 

dtS Plm (t) = 1 I dfi YC m [v ■ V(AnPo) + A Q P (V ■ v ± ) + r 2 (V ± (9 r P ) ■ Vx)v r ] I . 

Next we give explicit expressions for coefficients Spi m {t) for toroidal and spheroidal oscillations modes (see Unno et al. 



(1979)). Spheroidal separation of variables gives the following expressions for components of the oscillation velocity in spherical 
coordinates 

f 



U(r) Yi m (6, (j>), Ve=e- 1 ™ t V(r)dgY lm (9,<t > ), v^e 



V(r) -r-zd+YmiOii 
sma 



(B4) 



where U and V are radial and transversal velocity amplitude respectively. Substituting the velocity components (B4) into 
formula fed), we get the coefficients dtSpi> m > (t) for spheroidal oscillations 



dtSpi'm* it) = 



dtt Yy m , 



U (Y lm d r A n P + d e Y lm dgdrPo + -i^m d+frPb) + 
V sin J 

- (-1(1 + 1) Y lm AaPo + d B Y lm d e A Q P + -\-d^Y m d^A n P 
r V sin 



The oscillation velocity components for toroidal oscillations are 



v r = 0, vg = e W(r) -—d</,Yi m (0,cl>), v$ 



- iwt W{r)dgYi m (6,<j>) 



(B5) 



(B6) 



where W is transversal velocity amplitude. With the use of formula (B6), the coefficients in the expansion of dtSP for toroidal 
oscillations are 



dQ Yu 



W 



1 



r sin 6* 



(d^Ym OgAnPo - d e Y lm d^AnPo) 



(B7) 



APPENDIX C: EXPRESSIONS FOR GOLDREICH- JULIAN CHARGE DENSITY 
FOR MODES (1,1), (2,0), (2,1), (3,3) 

Pgj for < 6 < vr 

ii -;^t 3 \/f BW cos (0) sin ( 6 ') sin(<^)r, 2 



Pgj = e 



8 C 71- 2 r 3 



(CI) 



Pgj for < 6 < § 

3%/5B Wr» 2 (16 r - 3r» + 15 cos(40)r» - 12 cos(2 6>) (-4r + r„)) 



20 -iwt 

Pgj = -e 



256 C7T2 r 4 a/ 1 



sin(9)^ r» 



Pg°, for § < < 7T 

20 _ - !U t3V5BWr 4 2 (16r-3r« + 15cos(40)r«-12cos(20) (-4r + r,)) 
Pgj — e 

ain(fl) 2 r, 



256C7T2 r 4 \/l 



(02) 



(C3) 



Pgj for < 9 < vr 



21 — iutt V 6 

Pgj = e - 



- BW sin(6>) sin(0) r, 3 (21 r + 19 r, + 45 cos(2 0) r» 



3 7 

32 C7T 2 r 2 



(04) 



Pg 3 j for < 6 < \ 



33 _^ t 3V35BWsin(0) J sin(3 0) r.s ^ c „ 
Pgj = e 



256c7ri rT A /l 



(l5r 2 + 2rr» + 6 cos(40)r* 5 



6 cos(2 0)r» (-5r + r„) - 6\/2 cos(3 6») r - r» + cos(2 0) r, 
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(05) 
(06) 
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pfj for f < (9 < 7T 



„33 

Pgj = - 



,. tot 3V35Biyain(y Bin(3flr.* (l5r . + 2rr , +6 



(C7) 



6 cos(2 <9)r» (-5 r + r») + 6 cos(3 0) r»* a/2 r - r» + cos(2<9) r» 



(C8) 
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